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Q-i Abstract 

<D. 

We study generic waves without rotational symmetry in (2+1) - 
dimensional noncommutative scalar field theory. In the representation 
chosen, the radial coordinate is naturally rendered discrete. Nonlocal- 
ity along this coordinate, induced by noncommutativity, accounts for 
the angular dependence of the fields. The exact form of standing and 
propagating waves on such a discrete space is found in terms of finite 
series. A precise correspondence is established between the degree of 
nonlocality and the angular momentum of a field configuration. At 
small distance no classical singularities appear, even at the location 
of the sources. At large radius one recovers the usual commutative 
behaviour. 



1 Introduction and Overview 



Over the years, various attempts were made to introduce field theories which 
are nonlocal over a continuous or a discrete space. The hopes related to such 
a program were rarely fulfilled, often due to technical difficulties or structural 
issues. 

A relatively recent example of a nonlocal field theory (FT) without ob- 
vious inconsistencies and which can be handled analytically is given by non- 
commutative (NC) field theory pQ. One assumes that fields live over a space 
which is not commutative; for a plane one has [x, y] = i6 ^ 0. In the most 
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frequent approach to NC FT, nonlocality is simply attributed to the infinite 
number of derivatives in the star-product induced by the NC algebra. 

Using operatorial methods however, nonlocality can be quantified in a 
more precise way. Namely, the harmonic oscillator basis is used to naturally 
discretize the radial coordinate. Then, NC is traded for nonlocality along 
this discretized coordinate. Finally, precisely this nonlocality turns out to 
account for the angular (nonradial) dependence of the fields. Our goal is to 
understand free radial and nonradial, local and nonlocal, oscillation and prop- 
agation phenomena in such a framework, with discreteness and nonlocality 
explicit at the level of the degrees of freedom but manageable analytically. 

In the regime exactly opposite to the one we consider, namely when the 
kinetic term is neglected, interesting approximate radially symmetric solu- 
tions were obtained for scalar NC FT [2]. Subsequent perturbation theory 
in | confirmed [3]- [7] their existence but did not take the kinetic term on 
NC space fully into account. In spite of extensive work [2]- [IE], including 
exact results obtained for nonscalar theories [12]- [17], no complete account 
of free propagation on a discrete space seems available. In particular, one 
needs explicit solutions without radial symmetry to include configurations 
with nonzero angular momentum. 

Our contribution is the following. We provide new exact finite series ex- 
pressions for radially and nonradially symmetric oscillations and propagating 
waves, include sources and display the commutative limit. Along these tech- 
nical results two significant conclusions arise. First, no classical divergences 
appear at small distance; fields stay finite even at the location of the sources. 
Second, field configurations turn out to be local when they have zero angu- 
lar momentum and nonlocal whenever their angular momentum is nonzero. 
Actually the angular momentum and the (appropriately defined) degree of 
nonlocality of a field configuration can be identified up to a multiplicative 
constant. This last statement is proved through an extension of Noethers's 
theorem to operatorially valued field configurations. 

The paper is organized as follows. Section 2 introduces scalar NC FT, 
the radial basis which discretizes the field equations of motion and the simple 
procedure which trades NC for nonlocality. The nature of planar NC angular 
momentum is clarified in Section 3 through a Noether-type argument for 
operatorially valued-fields. A precise connection is established between the 
degree of nonlocality and the angular momentum of a field configuration. In 
Section 4 exact solutions are derived for the discrete field equations in terms 
of finite series. This is the main technical part of the paper. In Section 
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5 sources are included - without singularities appearing. The commutative 
limit is taken in Section 6, clarifying the nature of the NC solutions derived 
previously and confirming the precise way in which nonlocality is related to 
nonvanishing angular momentum. The Appendix collects known and less- 
known formulae required for the technical proofs. 

The general notation used is quite standard: n, k represent non-negative 
integers, the set of which is occasionally denoted by IV; the gamma function 
obeys T(z + 1) = zT(z); the binomial coefficient is C„ = , J^ !fc! ■ A discrete 
version of the logarithm, H n = 1 + 1/2 + .. . + 1/n, plays an important role 
in what follows. b~ n ^ = 5(n, k) is the Kronecker symbol, which is one when 
n = k and zero otherwise. For a local field 0, 4>(n) = n denotes its value at 
the point n. For a bilocal field $ that depends on two points, $(n, n') = Q n y 
denotes its value when the first point is taken at n and the second at n' . 



2 Equations of motion 

Consider a (2+l)-dimensional scalar field 0, depending on space coordinates 
forming a Heisenberg algebra (time is commutative and remains continuous): 

$(t,x,y), [x,y]=i61. (1) 

9 is a constant with the dimensionality of an area. The scalar field is 
consequently a time-dependent operator acting on the Hilbert space H. on 
which the algebra is represented. Since ([T]) implies 

[x, 4>(x, y)\ = i9^, [y, <f>(x, y)\ = -W^r> 

the field action, written in operatorial form, reads 

S = J dt Tr n (i ^ + l - [x, <ft\ % 0] + \ [x, f ] [y, 0] + V{^<t>)) (2) 

allowing for classical solutions <fi which are non-Hermitian, ^ ^ cf). We take 
V(</>V) = here. The equations of motion for the field are then 

i + l[x,[x,<P)) + ±[y,[y,<!>)} = 0. (3) 

In Cartesian coordinates, the solutions are plane waves 

~ e *(*-«+*»fi)-*-* j k 2 x + k 2 y = u 2 , (4) 
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formally identical to the commutative ones. In fact, due to the operators in 
the exponent, such waves have bilocal character [T9], in agreement with the 
considerations to follow. If a mass term m 2 <ft(f)/2 is inserted in Eq. fl2]), uj 2 
should be replaced by u 2 — m 2 . 

If the physical situation requires polar coordinates (a source emiting ra- 
diation, a circular membrane oscillating), then the harmonic oscillator, or 
radial, basis {|n)} 

N\n)—n\n), N = a^a, a = —={x + iy), n = 0,1,2,... (5) 

v 29 

is the natural one, cf. e.g. [10], and the equations of motion become 

0+^[a,[at,0]] = O. (6) 

N = K 5 -^ — 1) is basically (half) the radius square operator, in units of 6. 
Its eigenvalues n in (JSJ correspond to discrete points with radius growing like 
[n ~ Ijg) for large n. The NC plane is realized via (jSJ) as the semi-infinite 
discrete space of the points labeled by n. 
Origin of the recurrence relation 

We can sandwich any equation containing the operatorial field 4>(x, y, t) 
between (n'| and |n) states, eliminating NC in this way. The resulting field 
(n'\ (p{t) \n) = (f> n ' tn (t) is indeed a commutative object but also a bilocal one, 
as it generically depends on two different points, n and n' . Consider an 
harmonic time dependence <p n ',nif) — z wt< Pn'.n and use a \n) = y/n\n — 1), 
a* \n) = \/n + 1 |n + 1). Eq. implies that the equation of motion for 4>n',n 
is 

Vri + Wn + l</v+i,n+i + vn 7 v / n0 n '-i,n-i - (n' + 1 + n - \)4> n ',n = 0. (7) 

Above, A = |w 2 . Eq. (JTJ) is a recurrence relation of order two, describing 
the radial classical dynamics of a field which lives on a discrete space. The 
initial angular dependence of <p(x, y, t), due to its dependence on both x and 
y, is not lost. It is encoded in the two-index structure of 4> n ^ n) as will be 
seen explicitely in what follows. The operator is reconstructed from the 
c-numbers n / jn via 

<\>= E <f>n',n(t) i n \ ■ ( 8 ) 

The bilocality of the fields <p n ' >n thus implies the nonlocality of <j). 
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Radial symmetry 

If however the field depends only on the combination of x and y given 
by N, = 4>(N), we have radial symmetry. In this case is diagonal in the 
\n) basis, (n'| |n) = (fi n ',nO~n',n , and we have a local field. Defining n>n = n 
for simplicity, Eq. ([7]) implies 

(n + l)0 n+ i + n0 n -i + (A — 2n - l)0 n = 0, n = 0,1,2,... (9) 

again with A = u 2 6/2. The expectation value n = (n| |n) characterizes 
= Y^=o \ n ) 4>n (n\ uniformly at radius squared n. No angular dependence 
appears anymore. If a single value no is nonzero, then |n ) <fi no ( n o\ describes 
a field located at uq. 

A more suggestive form 

Returning to the general case (jSJ), consider the situation in which the first 
index of n ',n is greater than the second, n' > n. Define 0^ m - ) = n ',n) with 
m = n' — n > 0. The classical equation of motion of 0^ m ) follows from (JTj), 

Vn + m + 1 Vn+T^S+a + >/n + m N /n^ ) 1 + (A - 2n - m- l)0i m) = 0. (10) 

The index m stays constant throughout the above second-order difference 
equation in n. To include the situation in which the second index is greater 
than the first, introduce 4>^~ m ^ = <f> n ,n'i again with m = n' — n > 0. 0^ _m - ) is 
easily shown to obey the same equation as (jffl- Consequently, the indepen- 
dent solutions of Eq. ( TlOl are sufficient to characterize completely both cp^ 1 ' 
and 4>^~ m \ provided their boundary conditions are assigned independently. 
To summarize, if the notation 

Kun 2 = tfinlniL}' W, = \n X ~ U 2 \, U u U 2 G IN (11) 

is adopted, the expression §E§ turns into 

oo oo oo oo 

= E a m E \ n + m ) e "Vi m) ("I + E & ™ E \n) e iwi 0l" m) (n + m\ . (12) 

m=0 n=0 m=0 n=0 

Eq. (TlTJT) implies that configurations with different m can be freely superposed 
in (TT2"j) . as underlined by the insertion of the coefficients a m and b m , deter- 
mined solely through initial/boundary conditions. In contrast to \n) (n| 
which associates a value 0^ to the point n, \n + m)(j)^ n ' (n\ associates a 
value 0^ m ) to the two points n + m and n. Thus, m is a measure of the der- 
ealization of the field configuration it characterizes. Further statements can 
be made about the m-expansion ffT2"j) without solving the classical equation 
of motion ([TO]) , as we show next. 
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3 Angular Momentum versus Nonlocality 



Ref. [19] showed that the Cartesian coordinates plane wave excitations @ of 
a planar NC FT relate to one- dimensional dipoles described by a position x 
and an extension 5x = 6p y , with p y the linear momentum in the y direction. 
One cannot use y as a second independent variable due to [x,y] ^ 0; the 
conjugate variable p y turns out to be the natural substitute. 

A similar picture can be developed for radial coordinates. The bilocal 
quantity (fffl = n / jn can be described by two variables: a discrete radius 
squared given by n (by 2n + 1 actually) and an 'extension' m = \n' — n\. 
The analogy with the plane wave scenario suggests that m is related to the 
quantity conjugated to the angle, i.e. to the planar angular momentum. 

To prove this, we adapt Noether's theorem to the operatorial set-up of 
4>(x, y) and obtain the expression for angular momentum in the x — y plane, 
J z = J. First we identify the generator of rotations in the NC plane. Using 

" * * op 1 ry^ 

e iaR 0e -iaR = Q+ Q] _ _ [R ^ ^ Q]] _ ^ ^ Q]]] + . . . (13) 

and recalling N = \{^—^- 1) = a^a, we see that 

e taN xe~ taN = +xcosa + y sin a, e *aN~ e -iaN _ _xsina + y cos a. (14) 

Consequently, N generates rotations in the x — y plane. The variation of the 
field under an infinitesimal rotation is then 

5(f) = e iaff (f)(x,y)e- ia ^ -<j)(x,y) ~ ia[N,(f\, if a -»> 0. (15) 

Due to the trace over the Hilbert space appearing in the field action 
remains invariant under such unitary transformations. One can therefore 
proceed and adapt the usual way of thinking of the Noether theorem to this 
(still classical although) operatorial set-up. The conserved charge associated 
to the invariance under rotations, the angular momemtum, turns out to be 

J z = Trn iftla^a, <j>]. (16) 

The commutator appearing in (fTBT) is particularly simple in two instances, 

a)a,\n + m) (n\ — m \n + m) (n\ (17) 
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and 



a^a, \n) (n + m\ 



—m \n) (n + m\ . (18) 

These two equations already show that states 'delocalized' by an amount 
m are expected to carry m units of angular momentum. Using them and 
denoting the angular momentum of a field configuration by J z [0] we obtain 
for the expressions entering (TT2"j) 



J-. 



Lneiv 



\n + m) <f) { n m) (n\ 



n=0 



J-, 



-muj ^2 [ 

n=0 



-m)]2 



(19) 



(20) 



E l n )^~ m) (n + m| 

Dividing by the normalization factor iV+ = X^o^™**] 2 , respectively by 
= S^o[0n~ m ' ) ] 2 5 leaves us with the results (+mw), respectively (-mw). 
No matter what the precise values of 0^ and 0^ _m ' ) are, the index m de- 
termines the angular momentum, which is related exclusively to the degree 
of derealization of a field configuration. The fact that iV+ and iV~ may be 
formally infinite [similarly to the situation of commutative 2D radial waves, 
or even of harmonic ID waves] is irrelevant, since the normalization factors 
cancel out in the final result. If and N~ are finite or properly regularized 
[e.g. through an upper cut-off n max , or a multiplicative factor W 2 < 1 sup- 
pressing (fi^ at large n, to be removed only at the end of the calculations] 
we can further define 

' E \n + m) 0i m) (n\ , ( " m) = E \ n ) 0i~ m) (n + ™>\ 



hm) 



nelV 



and suggestively write the general solution ( 112"]) as an m-expansion 

= E K0 (+m) + M> ( ~ m) ]- 



(21) 



(22) 



meJV 



The operatorial solution 0( +m ) has derealization m, cf. ([21]) . and angular 
momentum +mco, cf. ( Fl9l) . Similarly, 0( _m ) has delocalization — m, cf. (121)1 . 
and angular momentum —moo, cf. ( 1201) . In consequence, (122 j) admits two 
equivalent interpretations. From the point of view of the theory defined on 
the discrete set of points n, it is an expansion in field configurations with 
well-defined nonlocality, more precisely bilocality, m. From the point of view 
of planar NC FT, Eq. ( 122]) takes into account the nonradial dependence of 
through an angular momentum expansion. 
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4 Bilocal waves via finite series 



In this section we solve Eq. ( TTUj) to obtain the exact form of (j)^ and 4>^~ m ^ . 
The difference equation (TTTJ|) describes travelling or standing waves on the 
semi-infinite discrete space of points n e IN. It is convenient to parametrize 
its two independent solutions as follows: 



(23) 



The functions /i(A) and /2(A), important for normalization, will be given 
later. They can be temporarely neglected since they do not depend on n. 

and which will turn out to be two polynomials in A, are denoted 
collectively by 0^ . In consequence, (123!) amounts for the time being to the 
substitution 



(m) _ /(n + m)!~ (m) 



nl 



(m) 



0. 



The field 0^ then satisfies the simpler recurrence 

(n + m + l)0i+\ + n$™\ + (A — 2n — m — 1) 
If the discrete derivative operator A is defined by 

A0 n = n+ l - 0„ 

and the shift operator E is defined by 

E<p n = n+ l, 

then the homogeneous difference equation (125]) can be rewritten as 



D 



(m) 



nA 2 E~ l + (m+ 1)A + A 



(m) 



0. 



annihilates the field 0^ m - ) , Vn G IV. 



The difference operator 
First solution 

Define for a real, the falling factorial power n~ , 



n 



T(n + 1) 
V(n+l-a) 



An" 



an 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 
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For natural a = k, n~ 



n[n 



(n — k + 1); this explains the name. 



We search for a solution of the form 

ft™) = a ™(a, \)n~ + <(a, \)n— + <(a, A)n 



j+2 



+ 



(30) 



i.e. an expansion in falling factorial powers of n. Equating to zero the 
coefficient of each n in 



D\ £ <(<t, A)n 

fc=0 



CT + fc 







(31) 



we obtain the indicial equation 

a(a + m) = (32) 

and the recurrence relation for the expansion coefficients a™ (a, A) 

{k + a){k + a + m) a r k n {a, A) + A a^a, A) = 0. (33) 

Eq. (1331) guarantees that ( 1301) is also an expansion in powers of A. For a — > 
we obtain the first finite series solution 



u 



(m) 



y> (~A) fc 
^ /c!(m + /c)! 



n! 



r(n + i} 



T(n+l-k) {n-k)\ 



n 



(34) 



From now on we drop the dimensionfull multiplicative constant ao by putting 
Oq = 1. It can be reintroduced whenever required. 
Second solution - infinite series 

The second solution cannot be obtained easily since we confront the case 
of roots differing by an integer in (|32"1) ; taking a = —m makes the coefficients 
diverge after some k. We therefore use the following procedure [20]: we solve 
the inhomogeneous equation 



k=0 



CT + fc 



a(a + m) n 



(35) 



take the derivative of its solution with respect to a, and then take the limit 
cr + m = e — > 0. This amounts to the evaluation of 



d_ 

fa 



E 



-\) k (a + m)T(l + a)Y{m + 1 + a)Y{n + 1 - k - a) 
,_ H T(k + l + a) T(k + m + 1 + a) T(n + 1 - k - a) 
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• (36) 



Applying d a to the summand in (1551) one obtains terms involving the digamma 
function = dl °^ x ^ = w^T- Using the poles-displaying expansion (IHTj) . 
then carefully observing the multiplicative cancellation between zeroes and 
poles in the limit o + m — > [and multiplying by a factor (— ) • {m — 1)! for 
convenience! we obtain the second solution 



w 



(in) 



E 



A m (-A)' 



Cn{H n _ k — H k — H k+m + 



m—l 



7) 



fe 4 (fc + m)! 

^ \ k ( m - k- 1)! ~ A fc+m (-)"n!(A;-n- 1)! 
^ k\(m-k + n)\ + k ^ +1 



(37) 



+ to)! 

As already mentioned, H k is a discrete version of the logarithmic function, 
# fe = l + ^ + i + ... + p fc = l,2,3... ; F = 0; (38) 

whereas 7 is the Euler-Mascheroni constant, 7 = limfc =00 (H k — lnk) ~ 0.5772. 
Second solution - finite series 

Although «4 m) solves ([25}, it is an infinite convergent series in A. To 
obtain a finite series, we search for a linear combination of wl m ) and ioi m \ 



4 m ) = a(A)4 m ) + 6(AVi m ), 
that is independent of A in n = and n = 1. If we choose to impose 



(m) , » , i in i ; , \ , 

t>o = a(A)«o + o(A)^o 



(m) 



,M 



we get 

a(A) 
KAl 



^! (m — k) ' m! Vm 



A r 

+ ^(- + 7 



0. 



00 yrn 

]^T{k + m)\ 



(39) 
(40) 

(41) 



Using (1371) . (1391 and (JUD to evaluate v[ m ^ we obtain, after a calculation 
featuring unexpected but essential simplifications, 

: (/(A)//, +0(A)(r, — -. (42) 



;i m) = a(A)4 m) +KAVi m) = ^r. 

m + 1 



Choosing 6(A) = e A ensures t> j m ' ) = ^W, which is constant (and convenient). 
In consequence, 

yk+m m \ 



^ X k ml 
. k =o k\(m-k) 



+ A r 



1 
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^ (k + m)\k 



u 



(in) 



(43) 
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We must still substitute in ( l4"5j) the expression (151]) for and the lenghty 
expression (1571) for u>^ m ). This leads to a long calculation involving products 
and sums of finite and infinite power series. Making use of the identities 
dZMZZD, (TTBIIBOD and dHSD we obtain the intermediate result 



,( m ) 



m— 1 



krik in— 1 



£A"(m-tf-l)!X; 



-A) fe C„ 



v=o 



fc=0 



(A; + m)! 



7V=0 



A? 



(m- JV-l)!(n + iV)! 



+ A^o(iV + m)! 



iV!(m + n)\ 



m QN 



+ 



(44) 



E^k+m f^lkl TT I TT TT TT \ _J_ \ " ^m+N-j 

fc=0 jtl ™ + J 



At this point, we are still left with a double sum in (144 p . Two cases appear: 
1. n < m — 1 . Then, the first term in (1441) can be written as 



m— 1 n 

E*V-"-i>.E (i+m)1 

n— 1 L m— 1 n ra+m— 1 n 

EE+EE + E E 

L=0 fc=0 L=n k=0 L=m k=L-m+l 



Wkr^k 

A) Vj n n < m —i 



(45) 



\L(_^k C k( m + k - L - 1 V- 



(m + k)\ 



2. n > m — 1 . Now, the same term in (1441) reads 



m— 1 



n / \ \krik 

j2\ N ( m -N-iy.J2 [ ] ' 



ym 

N=0 
rn—l L n— 1 



n>m— 1 



^ (A; + m)\ 



(46) 



L n+m— 1 n 

EE+E E + E E 

L=0 fc=0 L=m k=L—m+l L=n k=L—m+l 



\L(\k C k( m + k-L- 1)! 



(m + A;)! J 

Using either ( I4"5]) or ( 146]) - according to the case considered, grouping terms 
according to their power of A and observing some cancellations granted by 
(177] IHTl 188]) . we arrive both in case ( 143]) and in case ( 14B1) [hence irrespective 
of the relation between n and ml, at 



,M 



n-l 



' n— L 



E(-a) l E 



-l£<s+L 



L=0 



' (m + s)(m + s + 1) • • • (m + s + L) 



(47) 



This is our final result for v^. The first u^'s can be calculated easily. As 



expected, v, 



(m) 



for any natural m; then 



(m) 



(m + l)!/m! 



(m) _ (3 + m) - A 
;2 ~~ (m + 2)!/m! : 



(48) 
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(m) 
V3 



(m 2 + 6m + 11) - (2m + 8)A + A 2 



, , , , . VmGl. (49) 

(m + 3)!/m! v ' 

Concerning the functions f\ and /2 in (|23|) . the continuum limit behaviour 
of and requires, cf. Eqs. (T70|l and (1721, 

/i(A) = Af / 2 (A) = A"*. (50) 

Putting together (1231 . ijSTfl) . (I34l) and (T4T1) we obtain our final result for the 
two finite series that solve (fTOl) : 



aw 




n + m)! 



n! 



^2(m) 



(n + m)!^ 1 , 
A"" 1 - 2J( 



fc=0 
'n—L 



-A)' 



(fc + m) 



ft! 



C s n +L (m + s - 1)! 



(51) 



(52) 



l=o U=i (m + s + L)! J 

The solutions (I51II52P are finite sums and are linearly independent, since their 
Casoratian (the discrete analogue of the Wronskian) is nonvanishing: 



Dl 



A/n + m+ lA/n + l 



(53) 



,2(m) 



The general solution of f JTOj) is a linear combination = ci0^ m -* + C2<; 
with coefficients C\ and C2 determined from appropriate boundary conditions. 
Since obeys the same equation (flDjl . its general form will also be a 

superposition of the solutions ( JoTT) and ( 1521) but with coefficients determined 
from independently assigned boundary conditions, 0^ _m - ) = d x (^ m ^ +c 2 2 ^ m ^. 

Introducing now the expressions found for </>W and in ([12]) pro- 

duces the general solution of the difference equation of motion (fTOj) . As 
already proved in Section 3, m (respectively — m) characterizes simultane- 
ously the degree of nonlocality and the angular momentum of each bilocal 
configuration \n + m) (j)^ 1 ' {n\ (respectively \n) 0^~ m - ) (n + m|) in f[T2]) . The 
coefficients a m and b m are determined by the imposed initial/boundary con- 



ditions, for instance a Tl 



and b m = if only a field configuration of 



angular momentum ma; is excited. 

Connection to Laguerre polynomials 

It is worthwhile to mention a simple mathematical connection. If one 



writes 



(m) 



n 



(n+m) 



■L^" 1 ) [in contradistinction to (1241) ] one obtains 



+ 1)L<£\ + (n + m)LS + (A - m — 2n — 1)L™ = 



0. 



(54) 
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This is precisely the recurrence relation for the Laguerre polynomials, and 
it has two independent solutions expressible as polynomials in A. The first, 
£i(m) _ ( n +^) ! u M^ gives the well-known Laguerre polynomials ( l84"j) . The 
second finite series solution of (l5~4l) provides some kind of "Laguerre polyno- 
mials of the second kind" [which obey the Laguerre recurrence relation 
in n, (15"4"|) . but not the Laguerre differential equation in A] related to our 
through L 2 ^ = ( n+ "V ! t;( m ). Hence, we also obtained in this subsection a sec- 
ond finite series solution of the Laguerre recurrence relation (1541) . apparently 
unavailable in the literature. 

Radially symmetric solutions 

If one has radial symmetry, cf) = <p(N), the relevant difference equation 
simplifies to the local form (Q. Its first solution [18] follows immediately by 
taking m = in fl3] 



k=0 



To obtain the second solution one is confronted when m = with an indicial 
equation a 2 = that has roots which are equal, not differing by an integer. 
In spite of slight differences in proceeding due to that aspect, the resulting 
infinite series solution 

is still the specialization of ( 1371) to m = 0, namely w n = w^ ^. The last term 
in the right hand side of (|56l) was missed in [18]. Being negligible in the 
commutative limit, this term does not affect the large distance behaviour. 

Searching for a polynomial solution v n = a(X)u n + b(X)w n that obeys 
v o = and v i = 1 one obtains the m = specialization of ( 144")) . 

rt-l (_\\L L 

V n = E CnCL(^L- fc + fl-„-Jfc - 2# fc ). (57) 

L=0 k=0 

The form (|57|) . which appeared first in [11] without proof, can be considerably 
simplified. Namely, using (|85|) one obtains 

n-l (n-N ( _\s—lris+N ~] 
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which could have been obtained directly by taking m = in our general 
result (I47p . Further use of (177)) allows to reach an even simpler expression 



n-1 



(-A) 



k n 



E 



k (~ik 



(59) 



fc=0 



k\ 



k + j ' 



Both forms (158)) and (159)) are useful: the first is handier for some analytical 
manipulations whereas the second permits faster numerical evaluations. 

5 Sources 

Sources can be taken into account by introducing an inhomogeneous term j n 



Vn + m + Wn + + mVn$S + (A-2n-m-l)^ m) = j n . (60) 



The solution of Eq. ( 160)) takes into account an arbitrary distribution of 
sources j n and is consequently denoted by ^"^[jVt]. It can be obtained 
through the linear superposition of solutions ^^[j8 n ^ no ] which solve Eq. 
( 160)) with sources j n = jS njTl0 localized at an arbitrary but single point n . 

If no = one has the most interesting case - a source at the origin, 
jn = jS n ,o- If j = \frn\\~ m l 2 , the difference equation (160)) is now solved 
precisely by 0^ m - ) {+c<p l ^ m \ c arbitrary). Indeed, our hardly won 0^ m -* does 
solve the homogeneous equation ( flQl) everywhere except at n — 0, simply due 
to the fact that there the difference equation becomes first order and admits 
only one solution, which is </>^ m ). If a source j = \pm\\~ m l 2 is added at 
n = however, 0^ m -* is the particular solution of the resulting inhomogeneous 
equation. This is in line with the fact that v n enters the description of radially 
propagating waves, which require a source at the origin. If we take j arbitrary 
in j n = jS n} o, the solution will be (we write j = jo for later convenience) 



Consider then the case of a source at no > 1 (a ring-like source). We 
search for a solution of (1601 with source j n = jS n>no . A straightforward 
adaptation to the discrete case of the method of variation of constants (or of 
the method of Green functions) suggests to search for 



in (HOD 




\m/2 



(61) 



^ m) L?'<WJ = ci(n)^ + c 2 (n)^ m \ 



(62) 
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with ci and c 2 jumping only between n and n + 1, namely 

ci(n + 1) - ci(n) = rfi^no.n, c 2 (n + 1) - c 2 (n) = c/ 2 5 n0in . (63) 
Eqs. (IBTij) and (ITDjl imply then 

• ^2(m) ^l(m) 
d = J ^ _ ^ J (g^ 

1 V^o + m + Ia/^TT^D™ ' 2 V^o+m + lV^TT^' 

where is the Casoration of n (m ' ) and n (m ' ) at n = n , cf. Eq. ( 1531) . The 
solution for a source j5 n ,n of intensity j and location n is consequently 

^ [J^n.no] = Cx^M + erf™ ~ J0 n , no (fl^tf™ ~ ti^tf^) , (65) 

with the step function 9 m ^ n = 9(m — n) defined as 

fl if m - n > . . 

°^-\ o ifm-n<0. (66) 

Finally, denoting by $^[0] = ci0^ m - ) + c 2 0^ m - ) the general solution in the 
absence of sources, we obtain the general solution with an arbitrary source 
distribution j n as 

*i m) tin] = *™ [0] " £ JnAno (^V^ ~ 4?™^) ■ (67) 
no 

The sum is taken over all points n with nonzero sources, j no ^ 0. If a source 
jo appears at the origin one notes that, due to 0q = an d 0c/"^ — 0> 



the n = contribution in (j67|) reproduces (jblj) . The general solution (|6 
does not display singularities, even at the location of the sources. 



6 Commutative limit : types of waves 

It is useful to see first what happens to the difference equation ( flQj) in the 

commutative limit. As 9 — > 0, one has n ~ n' ~ |g — > oo and A = )■ 0, 

but A • n ~ (-f 1 ) 2 is finite; m stays finite as well. In this limit the difference 
operator applied in (flQj) to a function of n, 4>l™\ becomes the m-th order 
Bessel operator applied to a function of r, call it f( m \r). Indeed, taking 
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n 



— > oo, expanding the square roots in (JTDj) to order O(^) and replacing ^ 



by ^ one obtains 



Recalling that A = and passing via n = ^ from a function of n to a 
function of becomes as — >■ 0: 

j2 J?(m) 1 Jf (m) 2 

^ + ^ + ^-^)/ M (^0. (69) 

This is precisely the Bessel equation of order m for a function of independent 
variable cur; The solutions of f lTUj) should therefore reduce at large distances 
to linear combinations of the cylindrical functions of order m, ( l82i) and ( l83lh 
Our solutions are consistent with the above limit. Indeed, as n ~ ^ — > oo, 

A = and A • n ~ (<f ) 2 , Eqs. (EI} and ([82]) imply that a properly 

normalized i^" 1 -* becomes, as a function of r, the m-th order Bessel function 
J m (ojr), 

^±P^ «f0 J m (cr). (70) 

In this way we also establish the function fi in (1231) to be A m//2 , up to multi- 
plicative factors which go to 1 when A goes to 0. For instance, ortogonality 
of the (pn™ 1 in ( I5TT) . seen as functions of A G [0, oo), requires further multipli- 
cation of the LHS of ( 1701) by e~2. Not being of immediate relevance for our 
purposes, such factors are omitted. 

For the second solution things are less immediate. It is convenient to first 
consider w^. When n — > oo and A — > various terms disappear. First, the 
last one in fl37|) vanishes, 

^ X n+S 1 e A 

<T,1 77 < >"0. 71 

~~ fr{ (s + n)\n + s ~ n v ' 




g (-\) n \ s (s-l)\n\ 



!i (s + n)\(s + n + m)! 



Second, an ty^ m ) have the same commutative limit, due to the fact that 
the ratio | in (T4T1) vanishes as A — )■ 0. Using this observation and Eqs. ([3j 
and (TH3T) we obtain 



(n + m)! ^ m) _^ + j^)^^ _ 2 7 - log(^ 2 /2)]. (72) 
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This also establishes the function / 2 in ( 12 3 p to be \~ m / 2 '•. 

The correspondence between NC and usual waves can now be established 
through their behaviour at large distances. Given that the m-th order Bessel 
function J m {ur) describes usual radially standing waves (oscillations), Eq. 
( J70|) implies that its counterpart <p]^ describes radially standing NC waves 
(oscillations). Moreover, J m (ur) and Y m (ur) can carry angular momentum 
um or — um in usual planar field theory, in perfect agrement with the fact 
that 0^ m ** enters (fT2|) in combination with either \n + m) (n| or \n) (n + m\. 

On the other hand, the first Hankel function of order m 



describes waves which propagate outward radially and rotate angularly with 
frequency plus or minus um [unless waves with m and — m dependence are 
superposed, e.g. to render the angular part of the wave standing]. In con- 
sequence, the linear combination of <p]^ and (f 2 ^ 71 ^ which tends to H^ n (ur) 
as 9 —7- will describe a NC wave radially propagating outwards towards 
n = oo and carrying angular momentum +mu or —mu [unless two waves 
with opposite m are superposed]. This combination is easily found to be 



and displays angular momentum mu when it combines with |n + m) (n\ and 
angular momentum —mu when it combines with \n) (n + m\. Similarly, the 
easily found linear combination of (f)^ 171 ' and 0^ m - ) which, as 9 — > 0, tends to 
H^ur) = J m (ur) — iY m (ur) will describe "radially collapsing" NC waves. 

7 Appendix - Useful Formulae 

We collect here various formulae useful for the mathematical manipulations 
in the paper. To date, Eqs. ([76]), flTTJ), flTJ), flM}, (jSTJ, ([89]) could not be 
found in the literature accessed by the author. The form of these [possibly 
new] identities was suggested by the problem studied. Their proof requires 
only elementary methods and, occasionally, some effort. 

Recall first that H n = l + | + ^ + + Then we have the known 



H^(ur) = J m {ur) + iY m (ur) 



(73) 




(74) 



identity 



E 




H, 



(75) 



k 



n 



k=l 
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together with two useful generalizations of it, 

§" k (k + l)(kl2)...(k + m) =Hn+m ~ Hm ' (76) 

n / \fc— 1 (~<k+m n nm 

V I ^n+m l_ _ \ ~> '- y m+n-j /«~\ 

One can translate the denominator k in ( 1751) by a positive integer amount p, 
to find the identity (this time one can extend the sum to k = 0) 

Its generalization is, again for p integer and p > 1, 

y H^g m! = (p-l)!(n + m)! _ 

k^o (k + p)(k + p + 1) . . . (k + p + m) (p + n + m)\ 

The use of C*+™ = E]=o C m-i+j C n-j allows for a further extension: 

n ( \s-lr~<s+L n—1 n—j I \s— l/~ts 

f^ 1 {s + p)---{s + p + m) L - 1+ ' J f^ 1 {s + p)---{s + p + my 

after which one can use (179|) to complete the evaluation of the sum. 

The digamma function ^>(x) = dl °^ x ^ = j&Q, nas ^ ne following useful 
properties [21] (2; complex but not a negative integer, n positive integer) 

The series expansion of the Bessel functions [21] is required while taking the 
commutative limit of 

00 /n\2k+m 
■U*) = E M ' (82) 

^ fc!(fc + m)! 
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while the Neumann functions Y m [2TJ are also needed for w^, 

*s / \ n t / \ / , z , ^ (m-k-l)\ fz\ 2k ~ m 
7tY m (z) = 2J m (z)( 1 + log-)-J2- £j -(2J 

oo , x 2fc+m IT i IT 

-U~) k (?) uTu ' J? - (83) 



fc=0 



The Laguerre polynomials are given by (21 



rw M y r(n + m + l) (-A) fe 
nV ; Y r(Jfe + m + l)Jfe!(n-Jfe)!' 1 ] 



Finally, we write down the identities ensuring that is a polynomial 
of order n — 1 in A. Ifm = the following suffices, for N > n + 1: 



N-n 



Kcirtv^ + ^ - Ml) = g LL_1_^_|^ (85) 

For m > 1 however, we have two cases. If N <n, introduce for convenience 
the notation 

N m (jN 

E ^N^wPni^N-k + H n -k ~ H k - H k+m ) + ^ = A™ N . (86) 

k=o j=i n + 3 

Then 

(~) 4m _ / ^fc-i^fc (k — N — 1)! 

(iV+^)! n,7V ~ fc= V+i " (m + fc)! (87) 

and 

K,n = 0. (88) 
On the other hand, if TV > n + 1 the relevant identity is 

n m siN N-n /~m+m 



fc=0 
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